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Abstract. This paper reports a detection of the 30 mi- 
cron eniission feature from the C-rich Asymptotic Giant 
Branch (AGB) star IRAS 03313+6058 based on the ISO- 
SWS observation. Modeling of the spectral energy distri- 
bution shows that this emission starts at about 20 micron 
and possibly extends to the limit (45 micron) of the ob- 
servation. By assuming MgS as the carrier, the number 
ratio of Sulfur atom in MgS to Hydrogen atom in total, 
n(S)/n(II), is 3.0x10"^ from model fitting. A compari- 
son of this emission feature is made with other AGB and 
post-AGB objects. 

Key words: Stars: AGB and post-AGB - circumstel- 
lar matter - stars: chemically peculiar - Stars: individual: 
IRAS 03313-^6058 



1. Introduction 

The 30 micron emission was first observed in the bright 
AGB star, IRC-l-10216 (Low et al 1973). Observations of 
several carbon stars by Goebel & Moseley (1985) found 
that the emission also occurs in AFGL 3068. Recent ISO 
observations of some carbon stars turned up a few more 
candidates: AFGL 2256, AFGL 2155 and IRC-h40540 (Ya- 
mamura et al. 1998). This feature is seen not only in 
extreme C-rich AGB stars, but also in C-rich proto- 
planetary nebulae (PPNe; called also post-AGB objects) 
and planetary nebulae (PNe). Omont et al. (1995) ob- 
served five C-rich PPN objects with an unidentified emis- 
sion feature at 21 micron in their IRAS LRS spectra 
(Kwok et al. 1989) and detected the 30 micron feature 
in all of them. In the case of PNe a similar feature is 
observed for IC418 and NGC6752 (Forrest et al. 1981). 
A few suggestions for the carrier of this emission have 
been proposed. Because the feature has never been de- 
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tected in 0-rich objects, and is the broadest known fea- 
ture in the mid-infrared, solid species that form in the 
absence of oxides are taken into consideration. Omont et 
al. (1995) suggested iron atoms bound to PAH molecules 
as a possible emitter for the 30 micron feature. But Chan 
et al. (1997) raise doubt about this suggestion after the 
detection of this feature in direction of the Galactic cen- 
ter, where the existence of PAH molecules is question- 
able. A more acceptable candidate is solid magnesium sul- 
fide (MgS), first suggested by Goebel & Moseley (1985). 
The laboratory spectra of MgS samples showed very good 
agreement of band turn-on and cut-off wavelengths, as 
well as the overall band shape, with the observed fea- 
ture seen in AFGL 3068 and IRC-)-10216. A reasonable fit 
was achieved to the 30 micron feature of the PPN object 
IRAS 22272-f 5435 (Szczerba et al. 1997) by means of MgS 
grains with a distribution of shapes. In addition, MgS is 
one of the molecules which condensate at low temperature 
when no oxides are present. This is consistent with the de- 
tection of the 30 micron feature only in objects with cold 
dust shells. 

IRAS 03313-1-6058 was classified as a candidate ex- 
treme carbon star based on the similarity of its IRAS LRS 
spectrum to those of AFGL 3068 and IRC-l-10216 (Volk et 
al. 1992). It has no optical counterpart in the POSS plates 
(Jiang & Hu 1992). In the near infrared, the source was 
detected at K=15.6 mag and was not detected in J and H 
bands with upper limit of magnitude 17 and 16, respec- 
tively (Jiang et al. 1997). The color index between 12 and 
25 micron, based on the IRAS PSC catalogue, indicates a 
color temperature of about 250 K. Therefore, this is an ob- 
ject with a cold and optically thick circumstellar envelope. 
Though the IRAS LRS type of this object is 22, no clear 
silicate feature is seen (note that Kwok et al. 1997 clas- 
sified the LRS spectrum of this object as unusual). The 
detection of HON (1-0) fine (Omont et al. 1993) indicated 
the possibility of a C-rich nature. 0-rich maser lines such 
as OH (Le Squeren et al. 1992, Gait et al. 1989), H2O 
(Wouterloot et al. 1993) or SiO (Jiang et al. 1996) have 
not been detected at all and this is an indirect indication 
of a C-rich circumstellar envelope. The CO (2-1) line pro- 
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file suggests that this object is a late-type star rather than 
a young stellar object, and gives an expansion velocity of 
the shell of 13.9 km/s (Volk et al. 1993). 

2. Observation and data reduction 

The spectroscopic observation was carried out by using 
the SWS spectrometer (de Graauw et al. 1996) of the In- 
frared Space Observatory (ISO) satellite on 31 July 1997 
with the fastest scan speed covering full wavelength range 
from 2.3 /im to 45 fim (AOTOl, speed 1). The achieved 
resolution is about 200 ^ 300 with S/N higher than 100 
at wavelength longer than 5 /im. 

The original pipeline data were corrected for dark cur- 
rent, up-down scan difference, flat-field and flux calibra- 
tion by using the SWS Interactive Analysis (lA) at MPE 
Garching |^ The deglitching and averaging to equidistant 
spectral point in wavelength across the scans and detec- 
tors was done using the ISAP package. 

The ISOCAM imaging was performed on 31 July 1997 
with the CAM camera in the mode AOTOl with the fil- 
ter LW3, centered at 15 micron. The data were reduced 
by using CIA (a joint development by the ESA Astro- 
physics Division and the ISOCAM Consortium led by the 
ISOCAM PI, C. Cesarsky, Direction des Sciences de la 
Matiere, C.E.A., France; Cesarsky et al 1996). The ob- 
ject is still point-like at this angular resolution of 1.5 arc- 
sec/pixel. The flux through the filter LW3 is 59.03 Jy as 
measured by aperture photometry method, about twice 
of the 12 iim flux (30.87 Jy) given by the IRAS PSC 
catalogue. This difference could be caused by the strong 
emission in the mid-infrared range of the spectrum or by 
the object being variable. However, the flux (~ 39 Jy) 
at 15 micron of the SWS spectrum, which was taken at 
the same day as the ISOCAM image, does not match the 
photometric result from the ISOCAM image while it is 
in rough agreement with the 12 /im IRAS flux. By calcu- 
lating the IRAS fluxes at 12 and 25 micron from SWS, 
correction factors of 1.11 and 1.06 should be applied to 
the corresponding SWS bands, respectively, to agree with 
the photometric data from IRAS PSC. After correction, 
the spectrum is smooth and the shape is similar to its 
IRAS LRS spectrum. In the same time, the IRAS LRS 
spectrum should be multiplied by a factor of 1.27 to agree 
with the IRAS PSC data. These factors are within the 
calibration uncertainties of the ISO-SWS and the IRAS 
LRS. Such agreement may show that the flux calibration 
of ISO-SWS is reliable, and that the object is non-variable 
(the variability index from IRAS observation is 50, at the 
border between variable and nonvariable indices). There- 
fore, we suspect that the discrepancy between ISO-SWS 
and ISO-CAM results is related to the uncertainties in 

^ We acknowledge support from the ISO Spectrometer Data 
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Fig. 1. Observed and modelled spectral energy distribu- 
tion. The dot represents the near infrared and IRAS pho- 
tometric results and the thin solid line is the ISO-SWS 
spectrum starting from band IE. The flux at 15 micron 
measured from the ISO-CAM LW3 image is given by an 
open circle. The thick line is calculated from the model 
with the values of the parameters listed in Table 1 and the 
long-dashed line represents the radiation spectrum from 
the central star which radiates as a blackbody of temper- 
ature 2500 K. 

the flux calibration from the image at LW3 band, but the 
reason for this is not clear to us. Anyway this ISO-CAM 
result is shown in Fig. 1 as an open circle. 

3. Modeling 

3.1. Spectral Energy Distribution 

The overall spectral energy distribution of this object is 
observationally determined from the near infrared to the 
far infrared. The K band magnitude at 2.2 micron is 15.6 
(Jiang et al. 1997), i.e., the flux is 3.7x10-^ Jy. The fluxes 
at 12, 25, 60 and 100 /xm are 30.87, 43.44, 15.08 and 
6.40 Jy, respectively, from the IRAS photometric observa- 
tions ( quality index = 3 in all four IRAS bands). Combi- 
nation of these photometric results with the ISO SWS 01 
spectrum defines the observed spectral energy distribu- 
tion. In Fig. 1, the observational data are shown where 
the dots represent the photometric results and the thin 
solid line shows the ISO SWS 01 spectrum. The spectrum 
of ISO-SWS from 2.3 micron to 3.52 micron (band lA, 
IB and ID of ISO-SWS) is not shown because it is too 
noisy. 

The model we used to fit the observational data is 
described in detail by Szczerba et al. (1997). In brief, 
the frequency-dependent radiative transfer equations are 
solved under the assumption of spherically-symmetric ge- 
ometry simultaneously with the thermal balance equation 
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Table 1. Model parameters 



parameter 



value 



d 

Thinner 
dinner 

M 

^dyn 

Dust-to-gas ratio 

a- 

a+ 



2500 K 
8000 Lo 
4.25 kpc 
760 K 

5.8 X 10"cm 



7.0 
8.0 
1.6 
5.0 



X 
X 
X 
X 

0.005 iim 
0.25 iiiii 



10"' Mo/yr 
10* 
10" 



3 



p (power-law index of size dist.) 3.5 



for a dusty envelope. The radiation of the central star is 

assumed to be a blackbody. The mass loss rate is taken to 
be constant and the envelope is assumed to expand at the 
velocity derived from the CO (2-1) line observation. The 
dust opacity is represented by amorphous carbon grains of 
AC type (Bussoletti et al. 1987, Rouleau & Martin 1991). 

The appropriate values of important parameters for 
fitting the observational data arc listed in Tabic 1. The 
symbols have their usualmeanings, but more details can 
be found in Szczerba et al. (1997). Dynamical time tdyn 
means the time required for the matter to reach the outer 
radius of the envelope. 

The result from model calculation with the values 
listed in Table 1 is plotted in Fig. 1 by a heavy solid line, 
while the long-dashed line represents the radiation from 
the central star, which radiates as a blackbody of effec- 
tive temperature 2500 K. The luminosity and distance of 
the star depend on each other, as well as the dust-to-gas 
mass ratio and mass loss rate assumed. We adopted the 
value of 8000 Lq for the luminosity. The corresponding 
distance is 4.25 kpc. Since the outer radius Rout is 7.0 
X 10^'' cm, the angular diameter of the circumstellar enve- 
lope is predicted to be 22 arcseconds at this distance. This 
size is big enough to be resolved by the ISOCAM at the 
resolution mode of 1.5". However, the ISOCAM observa- 
tions centered at 15 /im are the most sensitive for the peak 
temperature of about 200 K which, according to modeling 
results, corresponds to radius of 4.1x10^^ cm and reflects 
a much smaller angular size of about 0.13 arcseconds. So 
the result of ISOCAM that no extension is seen may be 
attributed to a large distance of the object. The value of 
Rout is determined from the model-fitting to the observa- 
tional results and its choice affects mainly the flux inten- 
sity in the mid- and far-infrared. A reasoanable flt can 
be achieved with Hout values ranging from 5.0 x lO-'^^'cm 
to 9.0 X lO^^cm assuming a constant mass loss rate. For 
increasing mass loss rate (which means density distribu- 
tion steeper than r~^), the outer radius should be larger 
to compensate for the smaller far infrared emission of the 



outer enevelope layers, but range of the allowed changes in 
Rout is quite similar. Note, however, such density distribu- 
tion cannot be much different from that corresponding to 
the constant mass loss rate due to strong constraints from 
the 60 and 100 /zm flux densities unless we assume the 
dust optical properties at far infrared wavelengths have a 
less steep slope than in the case of amorphous carbon used 
here. 

The adopted dust-to-gas mass ratio of 5.0 x 10~^ cor- 
responds to mass loss rate of 8.0 x 10~^ Mq /yr. This mass 
loss rate is higher than 5.8 x 10^^ Mq/jt deduced from 
CO (1-0) line (Loup et al. 1993) or 1.4 x 10"^ Mg/yr in- 
ferred from CO (2-1) hue (Omont et al. 1993). Note, in ad- 
dition, that to get velocity of the shell aromid 14km/s dy- 
namical considerations would suggest even smaller value 
of dust-to-gas mass ratio (see Steffen et al. 1998), and 
in consequence a larger mass loss rate and a larger total 
mass of the envelope. However, the assumption of a den- 
sity distribution slightly more steep than r~^ would cancel 
the increase in total mass of the envelope. The calcula- 
tion of the mass loss rate from the CO line suffers mainly 
from the uncertainty of the distance and mass fraction 
factor for CO molecules. For example, the mass loss rate 
of W Hya derived from infrared water lines lies a factor 
about 30 above the estimates based on the CO line obser- 
vation (Neufeld et al. 1996). The case of Y CVn is similar 
in that the mass loss rate derived from interpretation of 
far infrared ISOPHOT images is 2 orders of magnitude 
higher than that foimd from the CO line (Izumiura et 
al. 1996). Izumiura et al.(1996) explained the result for 
Y CVn by suggesting that the far infrared and CO ob- 
servations represent different epoches of mass loss. There 
may be another possibility that, the mass fraction of CO 
molecules is overestimated so that the mass loss rate is 
underestimated. The mass loss rate of IRAS 03313+6058 
from modeling lies a little above that estimated from the 
flux at 60 micron, 4.8 x 10"^ Mg/yr (Omont et al. 1993). 
Derivation of mass loss rate from the flux at IRAS bands 
depends on the distance and on the bolometric correction 
factor which would induce some uncertainty. Though the 
mass loss rate derived from our modeling depends on the 
value of dust-to-gas ratio and density distribution, the 
result is relatively stable (probably better than a factor of 
two) because of the constraints required to fit the spectral 
energy distribution over the wide-range of wavelengths. 
The object then experiences quite strong wind, and has 
a circumstellar envelope perhaps as massive as 1.282 Mg. 
By considering that the outer radius Rout can vary in the 
range from 5.0 x 10^'^ to 9.0 x lO^^cm, the mass of the 
circumstellar envelope may be in the range of 0.92 Mq to 
1.65 M0 under the assumed dust-to-gas ratio of 0.005. 
Since the mass of the cirsumstellar envelope appears to 
be above one solar mass, the star could very possibly be 
an intermediate-mass AGB star. 
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3.2. The 30 micron emission 

As can be seen in Fig. 1, the object shows emission around 
30 micron which is superimposed on the continuum radi- 
ation from the star and the circumstellar envelope. The 
emission starts at 20 micron, peaks at about 30 micron 
and extends to longward of 40 micron. Because of another 
unidentified emission around 43.6 micron (see Discussion), 
it is difficult to define the cut off position of the 30 micron 
emission band from this spectrum, though the emission 
may extend to the long-wavelength limit of ISO-SWS. 

As described in the Introduction, MgS is regarded as 
a reasonable candidate to be the carrier of the 30 micron 
emission. We tried to model this with the optical constants 
taken from the tables based on laboratory measurements 
of aMgS(90%)-FeS(10%) mixture (Begemannet al. 1994). 
For our computations, we used two shapes for the grains, 
i.e. the CDE (Continuous Distribution of Ellipsoids) and 
Mie theory (spherical grains). Assumptions and method 
used are described in detail by Szczerba et al. (1997). 

Because the temperature structure of the circumstel- 
lar envelope is determined from modeling of the spectral 
energy distribution, and because there is little difference 
in the dust temperature structure between the largest and 
smallest AC grains (a_ = 0.005 and a-|- = 0.25 /xm) used in 
the calculations, the only free parameter after taking MgS 
into account is the number ratio between Sulfur atoms 
in MgS and total Hydrogen atoms n(S)/n(H). Under the 
CDE approximation a value of 3.0 x 10~^ for n(S)/n(H) 
gives a fit to the observed feature at the long wavelength 
wing, though there is a little inadequacy in the short- 
wavelength wing of the band. On the other hand, the Mie 
theory (spherical) grains can make up the emission at the 
short-wavelength wing. This means that a combination of 
MgS grains shapes may account for the observed emission. 
In Fig. 2, the observation and modeling of this band are 
shown. 



4. Discussion 
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Fig. 2. Fitting to the 30 micron emission. The thin solid 
line is the ISO-SWS spectrum and the thick line stands 

for the model fitting result to spectral energy distribution, 
the same as in Fig.l. The dot and dash lines represent the 
fitting by assuming the CDE and Mie approximation for 
the computation of optical properties of the MgS grains, 
respectively; n(S)/n(H) = 3.0 x 10"^. 
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It is interesting to compare the 30 micron feature 
of this object with those of other AGB and post- 
AGB stars. Fig. 3 exhibits the spectrum of IRAS 
03313+6058 together with the profiles of this feature 
for another AGB star, AFGL3068, and for the post- 
AGB object IRAS 22272+5435. The data for AFGL3068 
are taken from Yamamura et al. (1998) while the 
data for the IRAS 22272+5435 are from Omont et al. 
(1995). Because these two objects are much brighter 
than IRAS 03313+6058, their spectra are scaled down- 
ward to agree at about 20 /um with the spectrum of 
IRAS 03313+6058. Besides that IRAS 22272+5435 has an- 
other feature at 21 micron, the most evident difference is 
that the emission of IRAS 03313+6058 is much weaker 
than that of IRAS 22272+5435. On the other hand, 
AFGL 3068 exhibits a similar strength of the 30 fim band 



Fig. 3. Comparison of the 30 micron band in 
IRAS 03313+6058 (thin solid line) with the post-AGB 
object IRAS 22272+5435 (dotted line) and AGB star 
AFGL 3068 (dashed line). The modeling continuum spec- 
trum at this band to IRAS 03313+6058 is represented by 
the thick solid line which lies below the emission features. 

as does IRAS 03313+6058 which seems, however, to show 
more fine structures in this wavelength range. 

The 30 micron emission for IRAS 22272+5435 was pre- 
viously modeled by wing of MgS grains (Szczerba et al. 
1997). That fit resulted in an estimate of n(S)/n(H) = 
4 X 10~^ for the highest dust temperature distribution and 
1.6 X 10~^ for the lowest dust temperature distribution 
considered. The value of 3.0 x 10"*^ for IRAS 03313+6058 
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is very close to the case of highest dust temperature distri- 
bution of IRAS 22272+5435. Then there is probably little 
difference in the sulfur abundance between these two ob- 
jects and the strength of the 30 micron emission band may 
be influenced more strongly by other factors. The temper- 
ature may be one of the important factors through the way 
that lower temperature favors the excitation of this band. 
As the dust temperature of post- AGB envelopes is gener- 
ally lower than that of AGB ones, this emission is stronger 
in post- AGB stars. For example IRAS 22272+5425 has 
colder dust than IRAS 03313+6058 and much stronger 
emission at this band. From laboratory experiment, MgS 
is one of the molecules which condenses at low tempera- 
ture in the environment without oxides. Up to now, this 
band emission is detected in the C-rich AGB stars only 
with cold and high optical depth dust shells. This may 
indicate the existence of a critical temperature to form 
the 30 micron emission band carrier, and that such low 
temperatures together with approporiate chemical condi- 
tions are only possible in the extreme carbon stars. On 
the other hand, it is still not clear if carbon stars with 
smaller mass loss rates could form a carrier of this band. 
The existence of the 30 /Ltm emission in AFGL2155 and 
IRC+40540 (Yamamura et al. 1998), neither of them was 
classified as extreme C-stars by Volk et al. (1992), sug- 
gests that formation of the approporiate chemical material 
is much more common and that higher envelope temper- 
ature in other carbon stars (not extreme ones) does not 
allow us to detect this emission feature. But it could be 
as well that some special chemical reactions responsible 
for the carrier formation of the 30 /im band are efficient 
only when temperature is enough low and/or the chemical 
composition is approporiate. Unfortunately, since without 
optical spectra the determination of the chemical abbun- 
dances is rather impossible other methods of investigation 
should be elaborated, and especially better statistics cre- 
ated by the ISO data could help solve the problem of the 
30 iiui carrier formation. 

Besides this feature around 30 micron, some other fea- 
tures, e.g. absorptions around 7.5 and 14 micron, emis- 
sions around 41 and 43.5 /xm in the SWS spectrum of this 
object are not discussed here; they are currently under in- 
vestigation. We note only that the absorptions are prob- 
ably related to the C2II2 and/or HCN molecular bands, 
while emissions could be related to the crystaline silicates 
(especially as the enstatite mass absorption coeflScient 
matches these emissions well - see Jager et al. 1998). If 
crystalline silicate emissions are confirmed then it will al- 
low to deduce an evolutionary status of IRAS 03313+6058 
which could bring some more information on the exciting 
transition phase between oxygen and carbon rich parts of 
the AGB evolution. 
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